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SUMMARY 

Coefficients are derived for equations expressing the lateral force and pitching 

Inan^ translation and angular perturbations from a nom- 

lly centered rotating shaft with respect to a stationary seal. 

coefficients for the lowest order and first derivative terms emerge as being 
significant and are of approximately the same order of magnitude as the fundamental 

batlfrTn r T/"" Second derivative. shearpe^L- 

bation, and entrance coefficient ’variation effects are adjudged to be small. 

of the investigation delineated in this report defines the coeffi- 
cients or the equations: 


T = Ae + Be + ce + Da + Ed + Fci 

F = a£ + be + ce + da + ed + fd 


( 1 ) 


tb^r assumptions utilized, mathematical means employed, and conclusions derived 
thereby will be the objective of this report. 

presented^^^^'"*'^ additional terms upon a typical rotordynamlc system are 


INTRODUCTION 


leading to increasing pressures, speeds, and temperatures in jet 
g nes, igh performance pumps, and turbomachinery for rocket engine applications 
have been widely noted in the literature, e.g. , Rothe (ref. 1). 

The Increased Incidence of stability problems as a function of increased power 
author^ phenomenon of "subsynchronous whirl;’ has been noted by many 


ietv design features strongly contributing to this phenomenon in a wide var- 

IrJl uM] turbomachinery are annular shaft seals. These seals are in gen- 
1 utilized to separate regions of high and low gas or liquid pressure. In addition 
o eir strong effects as pseudo bearings and destabilizing devices, they offer, when 

*A11 work performed under NASA Marshall Space Flight Center Contract NAS8-27980 
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properly configured, in many designs the only real possibility for the introduction of 
damping in rotating systems. Various authors have delineated the extent of both the 
problems and remedies surrounding these possibilities; for example, Alford (ref. 2), 

Ek (ref. 3), Childs (ref. 4), and Gunter (ref. 5). 

A very significant part of the understanding of the dynamic behavior of annular 
smooth shaft seals has come from the pioneering work of Henry Black of Heriott Watt 
University, Edinburgh, United Kingdom, who in a series of papers (refs. 6-12) defined 
the effect of shaft displacements with both long and short seal assumptions up to and 
including the second derivative of the displacement. A number of workers have ex- 
tended his work significantly, these being Jensen (refs. 7 and 8), Hirs (ref. 13), 
Childs (ref. 4), and Alliare (refs. 5 and 14). The recent trend is to shift to a mod- 
ification of the bulk flow theory to define more closely a set of equations approxi- 
mating experimental results iHirs (ref, 13)]. 

On many significant recent problems, e.g., the problem on the high-speed rota- 
ting machinery for the Space Shuttle Main Engine (SSME) lEk (ref. 3)], Black’s equa- 
tions were the only analytical representa- 
tion available for significant computer 
modeling. Present efforts to extend this 
work include continuing work by Allaire 
(ref. 14) at the University of Virginia, 

Childs (ref. 4) at Texas A&M, and Fleming 
(ref. 12) of NASA Lewis, as well as exper- 
imental work with liquid oxygen and hydro- 
gen at the Rocketdyne Division of Rockwell 
International. 


This study represents a departure 
from the general approach being carried 
on by other investigators in that it ex- 
plores the effect of pitching moments and angulation in the seal, as shown schemati- 
cally in Fig. 1. 




FLOW 



Figure 1. - Orientation of references 
axes for seal equations* 


NOMENCLATURE 


The nomenclature in general is that employed by Black (ref. 6) and specific 
nomenclature is introduced as comments in the body of the nomenclature table, in the 
figures, or presented below. Additional functions are delineated in the text as 
required. 


A-F 

a-f 


} 


b 

D 


E 

F 


I 

K 


coefficients in force equations 

gap thickness, y(x,t), inches 

d/dt, differential operation 
with regard to time 

shaft elastic modulus 

force on seal, total due to 
perturbations 

shaft area moment of inertia, m ^ 
(in.^) 

effective bearing spring rate, 
N/m (Ib/in.) 


= force on seal, n is defined in the 
text 

= PR7TC/6X(l+5+2a) , general coeffi- 
cient for moment and force equa- 
tions 5 = 0.5 

= PRiraL/Xy 

= 3F/3c, N/m^ (Ib/in.^) 

= kinematic viscosity, m^/s 
(in. 2/sec) 
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L 

I = 

N 

P 

P 

R 

R 

a 

R 

r 

T 

t 

u = 

V 

W 

X = 


y,z = 

a 

Y 

6 = 
e 


In 

1 . 

2 . 

3. 

4. 

5. 


6 . 


7. 


length of seal, m (in.) X 

component spacing on rotor, m (in.) 

speed of rotation, rpm 

total pressure difference across 
seal. Pa (Ib/in.^) 

pressure within passage, p(x,t). 

Pa (lb/in.2) 


shaft radius, m (in.) 

axial Reynolds No. ZVy^/v = 

rotating Reynolds No. Ro)yo/v ^ 

passage time, L/V, seconds 

time, variable, seconds C 

axial fluid velocity, u(x,t), m/s 
(in. /sec) p 


mean fluid velocity, m/s (in. /sec) 

weight of rotor, kg (lb) ^ 

T “ 

axial distance, independent 
variable, m (in.) 

radial displacement of rotor mass ^ _ 

angular displacement ^ 

1.5 + 2a, ^ = 0.5 = 

weight density, kg/m^ (Ib/m^) X = 

lateral displacement, Y(t), m (in.) ^ 


( 1 ) 

0) 


friction coefficient, Ap = ^ Y— 


j. = area b 

n wetted perimeter, 2 

2 3/8 

a J 

loss factor per Yamada (ref. 15) 

functions of a, coefficients for 
the perturbation terms 

absolute viscosity, Pa-s (Ibf- 

sec/in.2 ) 

entry loss coefficient =0.5 
nominally 

mass density, -, kg-s2/m^ (Ib- 
sec^/in.'^) ^ 

XL/yo 

pitching moment taken about 
L 

X = -, N-m (lb-in.) 

pitching moment at x = 0, N-m 
(lb-in.) 

3F/3a, N/rad (Ib/rad) 

3t/3g, N (in.-lb/in.) 

3t/3a, N-m (in.-lb/rad) 
circular frequency, rad/sec 
design speed, rad/sec 


ASSUMPTIONS 

the following equations, assumptions are made as follows: 


The system (shaft within seal) is assumed to be centered. 

Perturbations about the nominal centered position are small. 

Whenever two second order terms are multiplied, they are assumed to be 
order, i.e., the product approaches zero. 


higher 


Short seal theory is assumed, i.e., perturbations have negligible effect in 
the tangential direction. 


The period of time required for a particle of fluid to travel from the inlet 
to the discharge is small compared to the frequency of the system. 

There is no change of properties in the working fluid as it passes through 
the process. 

The entrance coefficient into the seal is assumed to be 0.5. 
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8. The effect of fluid shear variation on the shaft moraent is neglected in the 
fundamental equations. 

9. Total linearity is assumed, i.e., any of the effects can be added individ- 
ually with no interdependence between individual effects. 

10. Tangential fluid velocity equal to one-half the shaft tangential velocity is 
developed immediately at the seal entrance. 

11. Perturbations are planar; however, Yamada's (ref. 15) work in defining X 
for rotating shaft effects applies. 

APPROACH 


The process of deriving the fundamental equations departs from that used by Black 
(ref. 6) in three important respects: (1) clearance is expressed as a function of 

both X and a, (2) the equations for continuity and velocity are completely rewritten 
as a consequence of (1), and (3) the order of integration is reversed. In addition, 
the first derivations lead to a result in which all perturbations, moments, and 
forces are taken at x = 0. This is done because of great simplification of the 
boundary conditions. These are then translated to perturbations about x = L/2, and in 
the case of the moments, is translated to t|l/ 2* Explanations of each step and 
necessary symbol derivations occur as needed in the text. 


DERIVATION OF FUNDAMENTAL EQUATIONS 


We first consider two plates of length 
L and unit width, nominally y in. apart, 
with the top plate fixed, the bottom plate 
perturbed: perturbations are positive up- 

ward and counterclockwise, as are forces 
and moments. Flow is positive to the right 
(Fig. 2). No fluid moves tangentially or 
normal to the surfaces of the plates. 

The steady-sta^ pressure distribu- 
tion through the fliw process is reflected 
by the following relationships: 



equations. 


p = Pi - ?2 = P(l-5 + 2a) |- 


( 2 ) 


P 


1 



(3) 


The fundamental equation for head loss through the flow process is 


9h _ X 

~ b g 


(4) 


From the Navier-Stokes 


equations for unsteady flow, 


_ — ^ ^ ^ 
p 9x ^ 9x 9t b 


( 5 ) 
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If 3. control volume of length dx, and bounded by the two plates, is now considered, 
the plates can be thought of as perturbed streamlines. The general continuity equa 
tion for an unsteady stream tube of fixed length ds is 


3(PA) ^ 8 


9 1 9s 

— ^ ^ 

Since p = constant, A=b, u A=ub, s=x 


+ - 5 — (pu • A) = 0 


or 


Now, approximately 


9b 9(ub) _ 

3t 3x ■ ^ 


, 3u 3b . 3 b ^ 


( 6 ) 


(7) 


( 8 ) 


or 


and U(t) = V + v(t) , 


9 u a . e . X . 
_ = — u + — + — a 

y^ y. 


o 


o o 


9u a _ e + xd 

9x y ^ ' y 

o ^o 


U = — (ex + ^ d) 

^o V 2 ; 


+ V + V + 


Vax 


(9) 


( 10 ) 


Higher order terms (HOT) are dropped. Now, integrate and set the following boundary 
conditions : 


and 


at X = 0, p = p^ = p^ - p(V^ + 2Va) 


at X = L, p = p. 


1 / . 2VA , x^ •• , / Vx , XVx^ \ . 

--P = — + — 2 e 

o n 'o y ' 


XV^x x^ 
+ -T- " ^ 6 y 


•• / Vx^ ^ XVx^ \ . 

- a + + ^ ) a 

V ^o 3y 2 / 

O 

) 


2 2 2 
V X . 3 XV X 


XV^x 


' ^o y 

o 


(11) 


Therefore 


1 S 2 

c(t) = +'~ 0 {v + 2Vv) 

^o 
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But 


= P = p(1.5 + 2a) f 


0 = V + - (1.5 + 2a)v + e + f 2(1 + a)e + \ 2oe 


La .. . V 


a + — 2(3 + a)a + — r 3(2 + 3a)a 


If now, p = Pgj.gg^y + Ap(t), then Eq. (11) becomes 




Vx ^ \Vk 


XV^x 

^o 


a + ^ 


2 3 

Vx ^ XVx 

3y„ 


/ V^x 

Uo 


2 2 

3 XV X 
2 2 
^o 


+ 1.5 Vv 


DERIVATION AND SUMMARY OF FORCE AND PITCHING 
MOMENT EQUATIONS 

The symbolism and form of the equations for a shaft is: 

T = Ae + Be + Ce + Da + Ed + Fa ) 

> 

F = ae + be + ce + da + ed + fa j 

The fundamental assumption of linear independence and superposition is made, i.e.; 

F^ = ae + be + ce \ 

F 2 = da + ed + fa I 

T, = Ae + Be + Ce V , , 

^ / (15) 

T 2 = Da + Ed + Fd I 

F - , T . + Tj ) 

In the above, all displacements, moments, and forces are taken about x = L/2. How- 
Gver, the first solutions will be taken about x = 0, and the results transferred to 

X = L/2. The subscript convention at x = 0 will be 


^o ^10 ^20 ’ ""o - "^10 '^20 
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For the cases about x = 0. certain simplifying parameters will be defined: the re- 

sulting form of the equations will be: 


20 


= d^a + e^a -f f^a = H (y a -b y Td -b y„ T^a) 


^10 ■ ‘^ o ^ ^ ^^ 10 ^ + 1^11 Te + y ^2 

2 .., 


21 


22 


■^10 " ^ o ^ ® o ^ ^ o ^’ " ^1 ^^^ 30 ^ ^31 + ^32 


(17) 

(18) 

(19) 

( 20 ) 


T^q - Da -b Ed -b Fd = (y^^a -b y^^ Td -b y^^ T^Ct) 

where T = L/V, and H, H^, and y^^j^ will be defined as part of the derivation. 

FORCE AND PITCHING MOMENT EQUATIONS FOR A SHAFT 

The definition of pressure perturbations seen in Eq. (13) is for two plates of 
unit width. Utilizing these, the following force and pitching equations can be writ- 
ten by integrating around a perturbed shaft as in Fig. 1. For example, 


F 
Rttp 


— = f Ap 

-np J ^ 


dx = 


j/v-b 2 VAv) 

1 ^ -b ^ e -b 1 

i VL^ 


A yo ) 

'2 6y ' 

0 

^2yo 

3^0 ) 


2 2 

+ + L 

O 2 ^ ^ 24y 
2y 


^ XVlS . 

^4yo nSyo 12yJ“ 


V^yo 2y^)“ 


-b 1.5 VvL 


( 21 ) 


RTrp 


2 


L . 

V 


2 r 2 

v-<-(2a + 3)v| + 1^1 ^ e -b ^ (3 -b 20)e -b 3ae 

I V 


H- 




24y^ 


2 a + ^ (8 -b 2a)d -b 12(1 -ba)aj 


and similarly. 


(vl!/ 

Rttp I 3 \ 


‘■V+20V+ (1.5)2v)+ 


) + !■'' fil-2.8L/l a\ 8o 
}* 8X [^2 ^ + V“(T + 4 )'= 




( 22 ) 


( 23 ) 


After eliminating v, v between these equations (23-24) and (12-13), the forces and 
torques can be defined. This is done piece-wise, then finally superposed. 
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FORCE DUE TO LATERAL DISPLACEMENTS 


This will be expressed in the form 


^10 = V + ^11 Te + Mi2 


(24) 


From Eq. (12), the equation expressing v,v as a function of e, t, e, is, with 
(a, d, a = 0) 


0 = |v + ^ (1-5 + 2a)v I + 


■"fx 


e + ^ 2(1 + a)e + ^ 2ae 


This can be expressed as 


0 = [TD + Y]v + ^ + 2(1 + a)TD + 2a]e 


L 2X 

Now, from Eq. (19), with a, d, d = 0 


10 

Rttp 


( 


2 

[TD + (2a + 3)]v + [T^D + (3 + 2o) TD + 3 ]e 


(25) 


(26) 


(27) 


Eliminating v, v between these two equations (26-27) with TD « Y, eliminating 
(TD) of third and higher power, 


’'in ■ » + U,, T4 + p„ T^fc-} 


10 


where 


10 


11 


22 


(28) 


oirRP 

6Xy 


, Y = (1-5 + 2a) 


10 


‘11 


y 


12 


= ^ 

Y 

= (9 + 12a + 4g^) y - 9a 

2 

Y 

- (8g^ + I8g^ + 18a) 


Y' 


This agrees with Black (ref. 6) except that in » Black has 19a; 18a is correct. 


FORCE DUE TO AN ANGULAR PERTURBATION 


This will be expressed in the form 


20 


= H (y^^a + y^^ Td + T^a) 


(29) 
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Proceeding as before, 


„ PR-rra , , 

" = ^ ’ Y = 1.5 + 2a 


_ 3L(2g + 6a + 3) 


22 


20 


1*21 = \ lY(2a^ + 11.5a + 12) - (6a^ + 18a + 9)] 
Y 


^ [Y (4.25 + 2a) - Y(2a + 11.5a + 12) + (6a^ + 18a + 9)] 


Y' 


PITCHING MOMENT DUE TO A TRANSLATION 

The pitching moment due to a translation will be expressed in the form 

''lO = "^o"^ ®o^ ^o'^ = ^1 ^»^30^ ^31 + ^32 


The equation expressing the boundary condition is (26) 

V o_ 

L 2A 


0 = [TD+yJv + 7 fr [T^D^ + 2(1 +a)TD + 2a]e 


From Eq. (21), a, d, a = 0, 


10 

RiTp 


VL 2 

(TD + 2a + 1.5 )v 


I 2 ^ 8 (i + f ) TD + |2j 


Eliminating v,v, as before. 


2 .., 


10 


■ “l <«30^ “31 ■“ “32 T 


where 


H, RTOLP , , 

1 = , Y = 1.5 + 2a 


30 2y 


1 f ^ + 3o; + 1 

4 2 


"31 Y I 3 

p =i 3\ _ 1_ (ol^ 3a 1\ 

^32 y\3 8/ 2V3 4 2/ 


2Y 

2 


Y 


2y' 


(30) 


(31) 


(32) 


(33) 
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The above pitching moment is about x = 0. 


TORQUE DUE TO AN ANGULAR PERTURBATION 

After proceeding as before, Eq. (20) is expressed in the form 

2 ,.. 


^20 ^ ^41 ^ ^42 ^ 


where 


„ prttol 

" "1^ 

2 

^40 " ^ 22 . 5 ] 


(34) 


= 2 I 39.75a + 33.75) - (22.5a + 50.67a + 22.5)] 

45y 

l2 2 2 

y ^2 = ~~3 1'^-25)y - Y(8a + 39.75a + 33.75) 

45y 

+ (22. 5a^ + 50.67a + 22.5)] 

TRANSLATION OF EQUATIONS TO x = L/2 


The translation of the solutions for perturbations about x = L/2 can now be 
written. These involve the a terms only. An x rotation at a = L/2, can be con- 
sidered as the sum of an a perturbation at x = 0 plus a displacement of -a L. , so 
that at L/2 there is a pure a perturbation. 2 


and 



(35) 

(36) 


These are both for a rotated at x = 1/2, but the moment is still at x = 0. The next 
step is to translate the resultant force and moment to x = L/2. The fundamental re- 
lationships are (F = F^) , 


T 




10 


+ T 


20 


(Fl t F^) 


(37) 


To summarize, these become: 

Forces ; 

F.,=aG + b£ + ce , since a = a,b = b,c = c 
1 o o o o o o 


F 


2 



(38) 
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Moments 




H 

O 

1 


- “c 1) " (»o - i)^ 


”^2 

T - 

20 

^2l=( 

”o ■ *o 1 - 1 + ‘‘o r)“ 



" i^o 

- ^1- 

o 2 o 4 / 



* t 

-ol- 

f ^ 

2 "■ "o Trr 

where 





a 

o 


, b 

0 

= H 

I , . H u^2 

d 

o 

and 

= «’^20 

, e 

o 

^ ^21 

T . f„ = H 

A 

o 

"l ^0 

, B 

o 

= »1^ 

bl ' '"o ° “l '^32 

D 

o 

= ^ ^0 

, E 

o 

= y 

‘41 • fo = "l ^2 


The symbols are previously defined. And finally. 


T = Ae + Be + ce + Da + Ed + Fa 
E = ae + be + ce + da + ed + fd 

where 


A = A - 


O % 2 ’ ^ ^ I , C = I 


D = D - A ^-d ^+a ^ 

o o 2 o 2 o 4 


o o 2 
2 


E = E - B y-e ^+b ^ 

o o 2 o 2 o 4 


L T T ^ 

F=F - C i-f ^+c 

o o 2 o 2 o 4 


and where 


a = a , b = b , c 


d = d - a ^ 
o o 2 


<1 JLJ 

, e = e - b - 
o o 2 


, f = - c 

o o 2 


c 

o 



( 39 ) 


( 40 ) 
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EVALUATION OF SEAL COEFFICIENTS 


The geometry and parameter values shown in Table I were used to evaluate a 
sample rotor. For comparative purposes, the properties of two fluids, water and 
steam under supercritical conditions were used. 

TABLE I. - GEOMETRY AND PARAMETERS 


SEAL 
R - 

RADIUS 

7.62 cm (3 INCHES) 

L - 

LENGTH 

2.54 cm (1 INCH) 

Vo - 

CLEARANCE 

0.254 mm (0.01 INCH) 

N - 

SHAFT SPEED 

30,000 RPM 

P - 

PRESSURE DROP 

344.6 N/cm^ (500 PSD 

^ - 

VISCOSITY (STEAM) 

2.067 X 10'^ N S/cm2 (3 X 10'^ LB-SEC/IN. 2) 

6 - 

DENSITY (STEAM) 

13.84 Kg/m^ (0.0005 LB/IN.^) 

^ - 

VISCOSITY (WATER) 

1.013 X 10'^ N-S/cm2 (1.47 X 10"^ LB-SEC/IN.^) 

6 - 

DENSITY (WATER) 

996.5 Kg/m^ (0.036 LB/IN.^) 

ROTOR 

E - ELASTIC MODULUS 

2.067 X 10^ N/cm2 (3 X 10^ LB/IN.^) 

1 - 

AREA MOMENT OF INERTIA 

24.97 cm^ (0.6 IN.^) 

K - 

BEARING SPRING RATE 

1.75 X 10® N/cm (10® LB/IN.) 

A - 

COMPONENT SPACING 

15.24 cm (6 IN.) 

w - 

ROTOR WEIGHT 

36.28 Kg (80 LB) 

G - 

i^K/EI 



As indicated in ref. 6, if all dynamic terms are retained in the coefficient 
formulation, the results are in the form of a frequency-dependent cubic divided by 
a first-order term. A Taylor's expansion was used for the denominator to obtain a 
simple polynomial form for the coefficients. Questions remain as to how many terms 
should be retained and what range of frequencies are allowed before the approxima- 
tion degrades. 

To answer the question of higher order terras, the equations for the seals 
were established in the cubic form and evaluated for the conditions of Table I. 

The resulting transfer functions are shown in Table 11, From these examples, it 
is obvious that the only affects that may legitimately be considered are those of 
stiffness and damping. Only those two affects were considered in the evaluation 
that follows. 

Seal coefficients were evaluated over a wide range of speed by assuming that 
pressure drop across the seal varied as the square of speed. The density for water 
was assumed to be constant but the density of steam was assumed to be proportional 
to pressure drop. The spring rates were found to vary with speed squared for both 
water and steam. The damping coefficients, however, varied with speed for the case 
of water and speed squared for steam. The restoring force and torque as a function 
of shaft translation and angulation are therefore described by 

/ \ - f ® U ^ I 
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table II. - SUMMARY OF SEAL COEFFICIENTS AT 30,000 rpm 


WATER 



9 IN/cm) = 7.84 X 10^ 

(’ ^ 548.3) (’ ^ 7M?) (’ i2l6o)/(’ * 7^) 

=5^7.84X10^+ 150.2S 

0 (N) = 6.806 X 10® 

(’ ^ ^ ’-948 * 45) 

=fc 6.606 X 10® + 264. IS 

X (Nl = -3.323 X 10^ 

(’ ^ 7^) iiiij/O * 7^) 

=:r -3.323 X 10^-31.43S 

0 (N-cm) « -1.2308X 10® 

(’ - 4^) (’ * sm) (’ ^ i4^)/(i * 7^) 

^-1.2308X 10®+ 19.76S 

STEAM 



6 (N/cm) = 9.327 X 10^ 

( 5335) (’ 48219) (’ * 85632)/(’ * ^iies) 

^9.327 X 10^+ 18.29S 

0 (N) = 1.0163 X 10® 

( 26295) (’ + _ 

;^^1.0163X 10® + 25.71S 

X (N) » -3.949 X 10^ 

(’ ^) (’ + «§«)/(’ 4^) 

^-3.949 X 10^ - 8.524S 

0 N-cm) = 3.536 X 10® 

( 41200) (’ 8O663) (’ * 233705)/(’ * 420w) 

^-3.536 X 10® + 2.688S 


where S * j(|) 



Where operating point variations are described by 


’ e 4)' 

-( 

\2 

a d 


if ' 

.X V'. 


03/ 

1 — 

> 

a 

1 

V' 

li)/ 


dt 

d_ 

dt 


d 1 

dt 

d 

dtJ 


(I) 


Steam 


■■b ^ 
dt 

B — 
L dt 




water 


ROTOR MODEL 


angulatlorat “anslation and 

the cantar. Saala wairplacjd « aLw^a"? a mass load at 

located halfway between each seal and ° shaft and an ideal bearing was 

model are contLnerirTabie I " the rotor 


The restoring force due to the shaft, 
lection (Y or Z) at the mass is given by 


bearings, and seals due to the lateral de- 


F/Y = N/D 

where 


(41) 


N = -2K 




+ 1 G_ 

2 iK 


(i> + X) 



(Q\l> 


- <t>x) 
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D.(l+f) + (f + |j0)|&+ (2+io) + 

12 '=) ^ + (3 + Ts)^®) W - w 

JL K 

Inserting values for the structural properties and for the frequency dependent journal 
coefficients, we obtain a dynamic force-deflection relation at the mass. 


F 

Y 


^ ^ -L 

^ ^ 72 

at 


d, + e, 4r + f 


1 1 dt 1 ^^2 


(42) 


Since the journal equations were developed for a condition where the flow is 
ostensibly axial, we assume that they are fixed to Couette coordinates. These 
equations must be transformed to rotor coordinates to determine the whirl orbit. 
Now assume, steady circular motion of the shaft in the y-z plane. Transforming the 
force-deflection equations to coordinates rotating with the shaft, we obtain 




0) 

2 ®1 


-4- 






(43) 


The 
the axis 


equations for the rotor with the central mass unbalanced a distance A 
of rotation in the y direction are 


from 


2 W A W 2 

- 0 ) — y = F + A — 0 ) 

g y g 


2 W 

- ( 1 ) — z = F 

g z 


(44) 


Solutions for shaft motion due to the unbalanced mass at a particular speed (o)) are 
obtained by first scaling the seal coefficients that were obtained at design speed 
(cu) . Then solutions of equations (41-44) produce the normalized orbit radius, or 
amplification ratio 


/ 2 . 2 
/ Y + Z 

A 

EVALUATION OF ROTOR RESPONSE 


Normalized rotor deflection vs shaft speed is shown in Fig. 3 and 4 for the 
steam and water journals, respectively. Rotor response considering only the trans- 
lation coefficient (a) is shown for reference. The rotor response with the complete 
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Figure 3. - Rotor response with 
steam as fluid. 


Figure 4. - Rotor response with 
water as fluid. 


outboard through the 

hv rh- ^ IS reversed. Reversed flow (inboard flow direction) was 

y aiiging tne sign of the off diagonal coefficients (cf) and ip) . 


seals and 
simulated 


a-ia '^“ evaluate sensitivity of the various coefficients, cases with the off- 
lagonal terms increased by 10 percent were run as well as a case where the dia- 
gonal coefficient (ip) was increased by 10 percent. A summary of these cases in 
terms of critical speed and effective damping (C/Ccr) are contained in Table III 


table III. - SUMMARY OF CRITICAL SPEED CALCULATIONS 
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DISCUSSION OF THE RESULTS 


When angulation affects are included, as well as the translation, the critical 
speed is shifted slightly and the whirl radius is significantly reduced. 

Off diagonal seal coefficients are very significant; without them, the effect 
of angulation is destabilizing. This indicates that the centers of pressure for 
rotation and for translation do not coincide. Reversing the journal flow (to the 
inboard direction) actually increased effective damping for the cases investigated. 

The second order "mass” affects of the journals were negligible and for whirl 
the mass affects may be deleted. 

For the two fluids that were considered in the journal, spring rates were very 
similar but damping coefficients were an order of magnitude higher with water. This 
is evidenced by an order of magnitude increase in damping at critical speed for tha1 
journal. For both fluids the spring rate varies with seal pressure drop and there- 
fore speed squared. The damping coefficient varies with speed for an incompressible 
fluid and the square of speed for a compressible fluid. 

CONCLUDING REMARKS 


From a study of the results, the following conclusions have been reached: 

1. Journal forces associated with the second time derivative of shaft motion 
should be ignored- The coefficients are very small and are significant 
only in a frequency range where the approximations required to obtain them 
are questionable. 

2. Forces associated with the first-time derivative are very important. 

3. Pitching moments and angulation effects are as significant as translation. 

4. The net effect of moment and angulation is a slight shift in critical speed 
and a significant decrease in peak amplitude. 

5. While the partial of moment with respect to angulation appears to be de- 
stabilizing, off diagonal coefficients result in net stabilization. This 
result is due to a lack of coincidence between the centers of pressure 
for translation and rotation. 

6. Seals, where the flow is outboard, showed lower damping than where the seal 
flow was inboard. Higher critical speeds, where a different mode shape is 
involved, may not exhibit the same trend. 

7. Damping forces at speeds below the design value are reduced more for a 
compressible fluid than for an incompressible one. 

8. Moments due to variation in fluid shear along the shaft surface are very 
small . 

9. Coefficient accuracy is not known in the absence of experimental data. 
Specifically this approach assumes full turbulent flow with a rotation at 
(jo/2 throughout the seal. Black (ref. 11) and Hirs (ref. 13) comment on 
this. For short seals, a significant error may occur because of the slow 
development of tangential velocity. 
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